Astronomy & Astrophysics manuscript no. aracil 


February 2, 2008 


(DOI: will be inserted by hand later) 





Metals in the Intergalactic Medium* 

Bastien Aracil^, Patrick Petitjean^'^, Christophe Pichon^, Jacqueline Bergeron^ 

^ Institut d'Astrophysique de Paris - CNRS, 98bis Boulevard Arago, F-75014 Paris, France, aracil@iap.fr 

LERMA, Observatoire de Paris-Meudon, 61 avenue de I'Observatoire, F-75014 Paris, France, petitjean@iap.fr 
^ Observatoire de Strasbourg, 11 rue de I'Universite, 67000 Strasbourg, France. 

Typeset February 2, 2008; Received / Accepted 

Abstract. We use high spectral resolution (R = 45000) and high signal-to-noise ratio (S/N ~ 35-70 per pixel) spectra of 19 
high-redshift (2.1 < Zcm < 3.2) quasars to investigate the metal content of the low-density intergalactic medium using pixel- 
by-pixel procedures. This high quality homogeneous survey gives the possibility to statistically search for metals at H I optical 
depths smaller than unity. We find that the gas is enriched in carbon and oxygen for neutral hydrogen optical depths thi > 1. 
Our observations strongly suggest that the C IV/H I ratio decreases with decreasing thi with log rciv ~ 1.3 x log thi — 3.2. 
We do not detect C IV absorption statistically associated with gas of thi < 1. However, we observe that a small fraction 
of the low density gas is associated with strong metal lines as a probable consequence of the IGM enrichment being highly 
inhomogeneous. We detect the presence of O VI down to thi ~ 0.2 with log tovi/thi ~ —2.0. We show that O VI absorption 
in the lowest density gas is located within ~300 km s^^ from strong H I lines. This suggests that this O VI phase may be part 
of winds flowing away from overdense regions. This effect is more important at the largest redshifts (z > 2.4). Therefore, at 
the limit of present surveys, the presence of metals in the underdense regions of the IGM is still to be demonstrated. 
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1. Introduction 

One of the key issues in observational cosmology is to un- 
derstand how and when star formation took place in the high 
redshift universe. In particular, it is not known when the first 
stars appeared or how they were spatially distributed. The di- 
i rect detection of these stars is challenging but the intergalactic 
■ medium (IGM) provides at least a record of stellar activity at 
these remote times. Indeed, metals are produced in stars and ex- 
pelled into the IGM by supernovae explosions and subsequent 
winds and/or by galaxy interactions. It is therefore crucial to 
observe the distribution of metals present in the IGM at high 
redshifts. 

The high-redshift intergalactic medium (IGM) is revealed by 
numerous H I absorption lines observed in the spectra of re- 
mote quasars (the so-called Lyman-a forest). It is believed that 
the gas in the IGM traces the potential wells of the dark mat- 
ter and its spatial structures: overdense sheets or filaments and 
underdense voids (e.g. Cen et al. 1994, Petitjean et al. 1995, 
Hernquist et al. 1996, Bi & Davidsen 1997). In the course 
of cosmic evolution, the gas is most likely metal enriched by 
winds flowing out from star-forming regions that are located 
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preferentially in the centre of massive halos. It is therefore not 
surprising to observe C IV absorption associated with most of 
the strong H I lines with log NiYi l) > 14.5 as these lines most 
likely trace filaments in which massive halos are embedded 
(Cowie et al. 1995, Tytler et al. 1995). The question of whether 
the gas filling the underdense space (the so-called voids) de- 
lineated by these overdense structures also contains metals or 
not is crucial. Indeed, it is improbable that winds from star- 
forming regions located in the filaments can pollute the voids 
entirely (Ferrara et al. 2000). Therefore, if metals are found in 
the gas filling the voids, then they must have been produced 
in the very early Universe by objects more of less uniformly 
spatially distributed. 

Absorptions arising through voids are mostly of low- 
column densities (typically of the order or less than 
A^(H l) = 10^^ cm^^). Given the expected low metalicities (typ- 
ically [C/H] < —2.5 relative to solar), direct detection of metals 
at such low neutral hydrogen optical depth is currently impos- 
sible due to the weakness of the expected metal absorption and 
statistical methods should be used instead. Lu et al. (1998) used 
the stacking method to increase the signal-to-noise ratio at the 
place where metal absorptions are expected and did not find any 
evidence for metals in the range lO^'^ < 7V(HI) < lO^"' cm~^. 
Although uncertainties in the position of the lines can lead to 
underestimate the absorption, they conclude that metalicity is 
smaller than 10^"^ solar in this gas. Note that this limit has 
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been confirmed by Ellison et al. (2000) using the same method. 
Cowie & Songaila (1998) introduced another method measur- 
ing the mean C IV optical depth corresponding to all pixels of 
the Lyman-a forest with similar H I optical depth (see Aguirre 
et al. 2002 for an extensive discussion of the method). They 
showed that the mean C IV optical depth correlates with thi for 
thi > 1. Ellison et al. (2000) used the same method on a spec- 
trum of very high signal to noise ratio and concluded that the 
data are consistent with an almost constant log C IV/H 1^—3 
down to THI ~ 2— 3. Using a large sample of high resolution 
data Schaye et al. (2003) find that the carbon abundance is 
spatially highly inhomogeneous. From simulations including 
strong assumptions on the UV background they conclude that 
the median metallicity is [C/H] = -SAltom- 
The method was applied to search for O VI by Songaila (1998) 
who found a significant amount of O VI in gas with a C IV 
optical depth greater than 0.05, and by Schaye et al. (2000) who 
claimed to have detected O VI in gas of mean H I optical depth 
as low as 0.1. However, Petitjean (2001) has shown that a non 
negligible fraction of the signal could originate in the wings 
of strong lines which are associated with overdense regions. In 
addition, note that three different groups have found different 
results concerning the nature of the O Vi phase (Carswell et al. 
2002, Simcoe et al. 2002, Bergeron et al. 2002). 

We have applied the methods introduced by Cowie & 
Songaila (1998), in a shghtly modified form, to a set of homo- 
geneous data of very high quality obtained in the course of the 
ESO Large Programme "The Cosmic Evolution of the IGM", 
in order to search for both C iv and O Vi absorptions in the 
IGM at 2 ~ 2.5. Section 2 describes the data and Section 3 the 
method, results of the investigation are presented in Section 4. 

2. Data 

The ESO-VLT Large Programme "The Cosmic Evolution of 
the IGM" has been devised to gather an homogeneous sample 
of lines of sight suitable for studying the Lyman-a forest in the 
redshift range 1.7—4.5. High spectral resolution (R ~ 45000), 
high signal-to-noise ratio (^35 and 70 per pixel at, respec- 
tively, 3500 and 6000 A) UVES spectra have been taken over 
the wavelength ranges 3 100-5400 and 5450-9000 A. Although 
the complete emission redshift range, 2.2—4.5, is covered, em- 
phasise is given to lower redshifts to take advantage of the very 
good sensitivity of UVES in the blue and of the fact that the 
Lyman-a forest is less blended and therefore easier to analyse. 
In particular, metal lines and amongst them the important O VI 
transitions can be more easily detected. 
Observations have been performed in service mode over 4 peri- 
ods (2 years). Details of data reduction and procedures used to 
normalise the spectra and preanalyse metal lines automatically 
will be described elsewhere. In brief, the data are reduced us- 
ing the UVES context of the ESO MIDAS data reduction pack- 
age applying the optimal extraction method and following the 
pipeline reduction step by step. The extraction sUt length is ad- 
justed to optimise the sky-background subtraction. The proce- 
dure systematically underestimates the sky-background signal 
but the final accuracy is better than 1%. Addition of individual 
exposures is performed using a sliding window and weighting 
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Fig. 1. Ratio of the true continuum to the recovered one used to correct 
the estimated continuum as a function of wavelength as measured in 
the rest-frame of the quasar. This is derived from normaUzation of 
simulated spectra. Different curves correspond to different redshifts 
(z = 2.0, 2.1, 2.3, 2.5, 2.7, 3.0, 3.2 from bottom to top). 



the signal by the total errors in each pixel. An automatic proce- 
dure estimates iteratively the continuum by minimising the sum 
of a regularisation term and a term, which is computed from 
the difference between the quasar spectrum and the continuum 
estimated during the previous iteration. Absorption lines are 
avoided using the estimated continuum. A few obvious defects 
are then corrected by hand adjusting the reference points of 
the fit. This happens to be important in small regions common 
to different observational settings and near the peak of strong 
emission lines. The automatic method works very well how- 
ever and manual intervention is at the minimum. We have care- 
fully calibrated this procedure using simulated quasar spectra 
(with emission and absorption lines) adding continuum modu- 
lations to mimic an imperfect correction of the blaze along the 
orders and noise to obtain a S/N ratio similar to that in the data. 
We noted that the procedure underestimates the true contin- 
uum in the Lyman-a forest by a quantity depending smoothly 
on the wavelength and the emission redshift and amounting to 
about 2% at z ~ 2.3. To calibrate this quantity we have simu- 
lated QSO absorption spectra drawing H I absorption lines at 
random from a population with the same column density and 
Doppler parameter distributions as observed, using the spectral 
resolution and noise as in the UVES data. We apply the same 
normalisation procedure and compute the ratio of the input to 
the output continua. This ratio is plotted in Fig. 1 for different 
emission redshifts. 

Here we use nineteen Unes of sight that are described in Table 1 
and Fig. 2. The mean redshift of the survey is approximately 
z ~ 2.4. 
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Table 1. List of lines of sight 



Name 




Forest 


Coverage 
C IV 


vr 


SNR 


HE1341 — 1020 


2.135 


1.55 


-2.10 


2.02—2.10 


2.00- 


-2.10 


66 


Q0122— 380 


2.190 


1.59 


-2.16 


2.07—2.16 


2.05- 


-2.16 


64 


PKS1448— 232 


2.220 


1.59 


-2.19 


2.07—2.19 


2.05- 


-2.19 


69 


PKS0237— 23 


2.222 


1.55 


-2.19 


2.02—2.19 


2.00- 


-2.19 


1 17 


HEOOO 1—2340 


2.263 


1.55 


-2.23 


2.02—2.23 


2.00- 


-2.23 


98 


Q0109— 3518 


2.404 


1.59 


-2.37 


2.07—2.37 


2.05- 


-2.37 


105 


Q1122 


2.410 


1.55 


-2.38 


2.02—2.38 


2.00- 


-2.38 


63 


HE2217— 2818 


2.414 


1.55 


-2.38 


2.02—2.38 


2.00- 


-2.38 


88 


Q0329 — 385 


2.435 


1.55 


-2.40 


l.Ql—l.AO 


2.00- 


-2.40 


61 


HE1158 — 1843 


2.449 


1.55 


-2.41 


2.02—2.41 


2.00- 


-2.41 


76 


n r, ] ,T+ / — Z.H-J 1 


2 611 


1.55 


-2.58 


9 09 9 


2.00- 


-2.58 


104 


Q0453-423 


2.658 


1.59 


-2.62 


2.07-2.62 


2.05- 


-2.62 


84 


PKS0329-255 


2.703 


1.62 


-2.67 


2.11-2.67 


2.09- 


-2.67 


45 


Q0002-422 


2.767 


1.63 


-2.73 


2.12-2.73 


2.10- 


-2.73 


69 


HE015 1-4326 


2.789 


1.63 


-2.75 


2.12-2.75 


2.11- 


-2.75 


97 


HE2347-4342 


2.871 


1.87 


-2.83 


2.40-2.83 


2.38- 


-2.83 


75 


HE0940-1050 


3.084 


1.96 


-3.04 


2.51-3.04 


2.49- 


-3.04 


84 


Q0420-388 


3.117 


2.09 


-3.08 


2.67-3.08 


2.64- 


-3.08 


88 


PKS2126-158 


3.280 


2.04 


-3.24 


2.60-3.24 


2.58- 


-3.24 


79 



The upper limit is the emission redshift blue-shifted by 3000 km s~ . The lower 
limits for the different redshift coverages {respectively, the forest, C IV and O Vi cover- 
ages) are the minimal observable redshift (due to instrumental limitation or the presence 
of a Lyman limit break) for Lymana, Lyman/? and O VI respectively. Lyman/3 is used for 
the C IV coverage to be sure that the H I optical depth associated to C iV is derived from 
at least two Lyman series lines. The SNR per pixel is averaged over the C IV regions. 
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Fig. 2. Limits of the Lyman-a forest along each line of sight in the 
survey. Different colors correspond to different redshift ranges. Light 
grey corresponds to the range between the Lyman-a and Lyman-/? 
emission lines where only Lyman-a absorption is expected. Dark grey 
corresponds to the range where Lyman-a is blended with other lines 
in the series. Black ticks correspond to the redshift at which Lyman/3 
is not redshifted enough to be observed. 



3. Pixel optical depth method 

In the highest SNR QSO spectrum obtained up to now, metals 
have been identified by direct detection of C IV absorption lines 
down to column densities as small as A^(C IV) = 10^^'^ cm^^ 
(Ellison et al. 2000). Most of these lines are associated with 
H I column densities larger than 10^^ cm~^ and are therefore 
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Fig. 3. The Lyman-Q forest (lower curve; here part of the 
HE 0151—4326 spectrum) is fitted consistently using the Lyman se- 
ries. The automatic procedure is therefore able to detect the presence 
of metals. The upper curve shows the fit of the forest where the proce- 
dure has recognised the presence of the strong C IV system indicated 
by tick marks. 



believed to trace overdense regions such as filaments of the 
IGM web. Our data have usually a detection limit of about 
10^^ cm^^ in the C IV region outside the Lyman-a forest. 
To test for the presence of metals in gas of smaller H I opti- 
cal depth, we use a variant of the pixel optical depth method 
(POD), that associates the H I optical depth with that in met- 
als (mainly C IV and O VI) on a pixel by pixel basis (Cowie & 
Songaila 1998, Ellison et al. 2000, Aguirre et al. 2002). 

The principle of the method is as follows. Once the Lyman-a 
forest is cleaned from all metals (see below), the H I Lyman-a 
optical depth is calculated in each pixel of the forest. When the 
Lyman-a line is saturated, other transition in the Lyman series 
are used. For each Lyman-a pixel, the observed optical depth in 
the metal transition at the same redshift is measured. Pixels in 
the Lyman-a forest are sorted according to their Lyman-a opti- 
cal depth and are gathered in optical depth bins that contains a 
large number of pixels. Then the median of the associated metal 
optical depths is calculated in the predefined bins. The median 
is chosen instead of the mean to avoid a few pixels with high 
optical depth to bias the mean optical depth. In that way it is 
possible to test for the overall presence of metals in low den- 
sity gas. The method has been extensively tested by Aguirre et 
al. (2002). Here we improve the method by carefully cleaning 
the Lyman-a forest for the presence of metals. 
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3. 1. Determination of mi and cleaning the Lyman-a 
forest from metals 

We want to estimate the Lyman-a optical depth in the largest 
number of pixels possible. For this we have to take into 
account that the absorption optical depth in any pixel can 
be due to several H i Lyman transitions at different red- 
shift and also to possible metal lines associated with strong 
systems. If Zem is the emission redshift of the quasar, and 
1215.67 and 1025.67 A are the wavelengths of the H i 
Lyman-a and Lyman-/? transitions, in the wavelength range 

[Amin,Amax] = [(l+Zem) X 1025.67,(l-HZem) X 1215.67], Only 

Lyman-a absorption is expected from H i. In the wavelength 
range A < (l+Zom)x 1025.67 the total absorption is due to the 
Lyman-a forest, the possible metal lines and also the absorp- 
tions in the different Lyman series transitions from the Lyman- 
a forest at higher redshift. 

We start considering A°J'pjjg^ in the wavelength range 

[Amin.Amax] = [(l+2:em) X 1025.67,(l+2;em) X 1215.67]. For 

each pixel, corresponding to A°['pj^,_^ = (1 + z)x 1215.67, we 
measure the effective optical depth in Lyman-a and calcu- 
late the corresponding effective optical depths at the positions, 
(1 + z) X Agj,j.;j,g, of all transitions in the Lyman series redshifted 
in the observed wavelength range, 

where /series is the transition oscillator strength and Ag(,j.igg the 
laboratory wavelength of the transition in the series. When the 
flux in a pixel is smaller than the noise rms at this position, we 
flag the pixel as a saturated pixel. The Lyman-a optical depth, 
Taipha in the considered pixel at Aj^j^pj^^^, is taken to be the small- 
est of the above Tgcrios that are larger than the noise rms at the 
corresponding position in the spectrum. This procedure allows 
us (i) to estimate the optical depth even if a transition is satu- 
rated; (ii) to avoid part of the blending effects. When all transi- 
tions are saturated, the pixel is flagged and is given a lower limit 
on T. Once Taipha is estimated over the above wavelength range 
[Amin, Amax], it is possiblc to subtract to the spectrum the corre- 
sponding optical depths in all Lyman transitions. Then it is pos- 
sible to go on with the determination of Taipha over the wave- 
length range [Ami„,Amax] = [(l+2em)x 1025.672/1215.67, 
(iH-Zem)x 1025.67] usiug the same procedure. 
The advantage of this procedure is not only to increase the red- 
shift range over which the study can be performed but to clean 
the Lyman-a forest from most of the polluting strong metal 
absorptions. Indeed, the strong metal absorptions are left over 
once the H i absorptions have been subtracted. Figure 3 illus- 
trates how successful the procedure can be in recognising auto- 
matically strong absorptions due to metals in the forest on the 
basis of the forest properties only. 

3.2. Cleaning the spectra 

The spectra have been scrutinised to flag all portions of the 
lines of sight where strong metal lines hide the information 
we need. These regions and their complements (O vi and H i 
wavelength ranges in the forest, C iv wavelength range in the 
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Fig. 4. Portion of the spectrum of HE 0151—4326 centered at 
0=2.4178 showing a strong Lyman-a absorption line with no C IV 
absorption associated but two satellites H I absorptions have strong as- 
sociated C IV and O vi absorptions. These satellites are narrow, there 
is no sign of strong displacement between H I, C iv and O VI, N V is 
absent and covering factor is unity: these lines are kept in the analysis 
and are not considered as "associated" with the quasar. 



red) are removed from the analysis. This includes in particu- 
lar strong Fe II and Mg II systems and all associated systems. 
The latter are characterised by partial covering factors and large 
NiC iv)/iV(H I), NiN v)/Ar(H i) and iV(0 vi)/Af(H i) column 
density ratios. When in doubt, we have kept the absorptions 
as in the system shown on figure 4. The latter is observed at 
Zabs = 2.4178 towards HE015 1-4326, 42600 km.s"! from 
the emission redshift (Zem = 2.789). 

This means that the C IV and O VI optical depths could be 
shghtly overestimated although we are confident that enough 
care has been exercised and no obvious system has been 
missed. We have also restricted the analysis to absorptions lo- 
cated at least 3000 km s~^ away from the emission redshift. 



3.3. Metal optical depth 

Once the optical depth in the Lyman-a forest, thi, is known, 
the corresponding optical depth of associated metal transitions 
at the same redshift can be derived. 

To avoid the effects of blending, we use doublets such as 
C ivAA1548,1550 and O viAA1031,1037 to secure the opti- 
cal depth. For each of the doublets, the first line (#1 of rest 
wavelength A?) has an oscillator strength twice larger than the 
second line (#2 of rest wavelength A2). Therefore, the true op- 
tical depth in transition #1 must be twice larger than the true 
optical depth in transition #2. We therefore can impose that the 
two optical depths be consistent. Using the observed optical 
depths at the position of the two transitions, rf^^ and T2^'^, we 
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Table 2. 


Determination of 


Tl — Tmctal 










Saturated 


Line 1 
In the noise 


OK 


Line 2 


Saturated 
In tlie Noise 
OK 


Lower limit 

tS""' X 




^obs 

Check 
Check 



want to determine the best estimator of the optical depth, ti, in 
the first transition. Starting from the smallest wavelength, we 
calculate for each pixel of wavelength Ai, the observed optical 
depth at this position t^^^ = — In(Fi) and the observed optical 
depth T2^^ = -ln(i^2) at a position A2 = Ai x A2/A". 
Two cases require special consideration. The potential absorp- 
tion is said to be "In the noise" if F > 1 — k, x a with a the 
rms of the noise in the corresponding pixel and k = 3. When 
F < K X a, the absorption is said to be "Saturated". 
We use for ri, either T°^' or t"^"" X 

fi I J 2 depending on the 
status of the absorptions. The different cases are summarised 
in Table 2. In the case the two lines are neither in the noise 
nor saturated (the lines are "OK"), or if #1 is "OIC' and #2 
is "In the noise" (cases denoted "Check" in the Table), some 
additional test has to be performed. First Fi is used to predict 
F2^^ that is the expected flux F2 in case the absorptions are 
only due to unblended lines from the doublet considered. If 
F2""P > (F2 - as), then n = T°b^ If not, n = t"""' x /1//2. 
Once Ti is defined and if Fi < 1 — k x cti, we replace T2^^ 
at position A2, by T2^^ — ti x /2/./1. This coiTesponds to a 
self-contamination correction. 

3.4. Simulations 

To check our two procedures, normalisation of the spectrum 
and recovery of the C IV optical depth, we simulated QSO 
spectra drawing H I absorption lines at random from a pop- 
ulation with the same column density and Doppler parameter 
distributions as observed, using the same spectral resolution as 
in the UVES data and adding noise. We use S/N = 10000 to 
simulate perfect data and S/N = 70 to simulate observations. 
We include C IV absorption at the same redshift as H I, assum- 
ing a fixed log r(C lv)/r(H l) ratio. We recover r(C IV) versus 
r(H I) and plot results in Fig. 5 for log r(C IV)/t(H l) = -2.5 
(squares), —3 (crosses), —3.5 (circles) and —10 (diamond) to 
simulate data without any enrichment. We applied the same 
pixel-to-pixel procedure after normalisation of the data either 
with the known continuum (top panels) or using our normal- 
isation procedure (bottom panels). We plot r(C IV) recovered 
by the procedure versus t(H l). It can be seen that we can re- 
cover coiTectly r(C IV) in all cases down to log thi ^ 0. The 
presence of metals is revealed by a correlation between thi and 
Tmctais and the plateau at small H I optical depths is a conse- 
quence of limited S/N ratio. It is clear that these simulations 
give us full confidence in our results. 

4. Results 

We apply the following procedure either for one line of sight or 
for the complete data set, combining all lines of sight. Pixels in 
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Fig. 5. Test of the normalisation and pixel-by-pixel procedures. We 
simulate QSO spectra drawing H I absorption lines at random from a 
population with the same column density and Doppler parameter dis- 
tributions as observed, using the same spectral resolution as the UVES 
data and adding noise, S/N = 10000 in the left panels and S/N = 70 in 
the right panels. We include C IV absorption at the same redshift as H I 
assuming some fixed ratio, log r(C IV)/r(H I) = —2.5 (squares), —3 
(crosses), —3.5 (triangles) and —10 (diamonds). We apply the same 
pixel-to-pixel procedure (all panels) to the data without normalisa- 
tion (the continuum is assumed to be known; upper panels) and with 
normalisation (using our automatic procedure) of the spectra (bottom 
panels). We plot r(C IV) recovered by the procedure versus r(H I). It 
can be seen that we can recover correctly r(C IV) in all cases down to 
log T ~ 0. 

the Lyman-a forest are sorted in order of increasing H I Lyman- 
a optical depth and are gathered in predefined bins. In each of 
the bins, the median of the corresponding Tmctai, measured at 
the same redshift, is then calculated. Using the median instead 
of the mean avoids the measurement to be biased by the pres- 
ence of a few strong absorptions. Individual plots for the nine- 
teen lines of sight are shown for C IV and O VI in Fig. 6 and 7, 
respectively. 

4. 1. Presence of C iv 

We plot in Fig. 8 the median of the C IV optical depth ver- 
sus the median of the H I optical depth when combining the 
whole sample. It can be seen that there is an excess of C IV 
for Thi > 1 . This excess is present for nearly all the individual 
lines of sight (see Fig. 6). For thi < 1, there is no detection of 
C IV. In particular, the apparent optical depth is statistically the 
same and consistent with zero for all the bins corresponding to 

log THI < 0. 

The thick dashed line corresponds to a fit to the data for 
log Thi > 0: log tqw = 1.3xlog thi — 3.2. The dotted 
diagonale line is drawn for illustration and corresponds to 
log(rciv/''Hi) = —2.6. It can be seen that the diagonale line 
is rejected by our data. Indeed, the present data are inconsistent 
with a uniform distribution of C IV in the Lyman-a forest with 
constant tqiv I ^hi ratio. On the contrary our results strongly 
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Fig. 6. C IV optical depth versus H I optical depth for individual quasars in the sample. 
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Fig. 7. O VI optical depth versus H I optical depth for the 19 individual quasars in the sample. 



suggest that the C IV/H I ratio increases with increasing thi- 
It is also apparent that for log thi ~ 0, log(C IV/H l) < -3.2. 
Although the results by Schaye et al. (2003) favor a slope closer 
to the diagonale than what we observe, their error bars are 
larger and their measurements are consistent with our findings. 

It is important to note also that the scatter in the C IV optical 
depth is very similar at different H I optical depths. This means 
that although the mean C iv optical depth decreases with de- 



creasing Thi, large C IV optical depths can be seen at nearly any 
H I optical depth. This is illustrated in Fig. 4 where a strong 
Lyman-a system is seen with no C IV absorption associated 
down to our detection limit ~10^^ cm^'^, but two satellite H I 
absorptions have strong associated C IV and O VI absorption 
(see also Bergeron et al. 2002). It must be noted that these sys- 
tems have quiet dynamics (the different absorptions of C IV, 
O VI and H I are fairly well centered on top of each others), they 
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Fig. 8. Median of C iv optical depth versus median of H I optical depth 
is plotted with black cross. Only pixels with velocity difference larger 
than 3000 km s~^ compared to the emission redshift of the quasar are 
used. Grey circles correspond to the previous finding by Ellison et al. 
(2000). The dotted line is drawn for illustration and corresponds to 
log C iv/H I ~ —2.6 . The dashed thick line is a fit to the data for 
log thi > 0; log rciv = 1.3 X log thi — 3.2. The horizontal dashed 
line corresponds to the median of rciv for log thi < — 1.3. 

do not show evidence for partial covering factors and have very 
little or no N V associated absorption. Therefore, their proper- 
ties strongly differ from those of systems associated with the 
quasar and they must be truly intervening systems. 

4.2. Presence of OVI 

Schaye et al. (2000) performed a similar pixel-by-pixel search 
for O VI absorption in eight high-quality quasar spectra span- 
ning the redshift range z = 2.0^.5. In the redshift range 
2 < z < 3, they stated to detect O vi in the form of a pos- 
itive correlation between the H l Lyman-a optical depth and 
the optical depth in the corresponding O VI pixel, down to 
Thi ~ 10~^. On the contrary, they did not detect O Vi at ^; > 3 
and considered this is consistent with the enhanced photoion- 
ization from a hardening of the UV background below 2; < 3, 
although this could also be caused by the high level of contam- 
ination from Lyman series lines. 

It can be seen from Fig. 9 that our data confirm the detection of 
O VI down to Thi ~ 0.2 that is at smaller H I optical depth than 
for C IV which is seen down to thi ~ 1 only. To investigate if 
the signal for 0.2 < thi < 1 comes from all parts of the spec- 
trum however, we have divided the sample into two subsam- 
ples. The first subsample contains all H i pixels located within 
Av from a strong absorption with thi > 4 and the second sub- 
sample contains all other pixels. We then vary Av. It can be 
seen on Fig. 10 that for Av < 300 km s^^, the signal signifi- 
cantly increases for the first subsample. This means that O vi is 
predominantly seen in the vicinity of strong Lyman-a absorp- 
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Fig. 9. The median O VI optical depth is plotted versus the median 
H I optical depth. The crosses are the observed data, other curves are 
simulations where the ratio log O Vi/H I is given fixed values, —5.0 
(no metals), —2.5, —2.0, —1.7, —1.5 from the curve labelled (a) to the 
one labelled (e), respectively 




-1.5- 



-1.0 -6.5 o!o ols 

log(THi) 

Fig. 10. Median O vi optical depth versus median H i for pixels lo- 
cated at less than 500 km s~^ (upper dashed curve) or less than 
300 km s~^ (upper dotted curve) from strong absorptions (thi > 4) 
and for pixels located at more than 500 km s~^ (lower dashed curve) 
or more than 300 km s~^ (lower dotted curve) from strong absorp- 
tions. 

tion lines. Indeed, the signal disappears for the second subsam- 
ple; there is no O vi absorption for 0.2 < thi < 1. This new 
result strongly suggests that the O vi absorption arises mostly 
in regions close in space to those responsible for strong H I ab- 
sorptions and a most likely explanation is that the O vi phase 
is part of winds flowing away from overdense regions. 
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Fig. 11. The C IV (top panel) and O VI (bottom panel) optical depth is 
plotted versus the H I optical depth for z < 2.4 (filled circles and solid 
line) and z > 2.4 (triangles and dashed line). 



To derive a quantitative limit on the amount of O VI present in 
the IGM, we have performed simulations of the metal enrich- 
ment of the Lyman-a forest. Artificial spectra are created draw- 
ing absorption lines at random from a population with column 
density and Doppler parameter distributions consistent with 
those observed (Petitjean et al. 1993, Hu et al. 1995, Kirkman 
and Tytler 1997, Kim et al. 2000). The number of lines is ad- 
justed so that the mean absorption of the simulated spectra is 
the same as the observed one. Although this is probably a rough 
assumption, constant O Vl/H I ratio is assumed and O VI ab- 
sorption is added accordingly. Noise consistent with the data is 
added and the whole procedure described in Section 3 is ap- 
plied to the simulated spectra. Results are plotted in Fig. 9. It 
can be seen that data are consistent with log O Vl/H I ~ —2, 
a value within the range found for individual O VI systems 
(Bergeron et al. 2002). 

4.3. Evolution with redsliift 

We have divided the sample in two redshift bins with similar 
number of pixels in each bin, z > and < 2.4. We plot in Fig. 1 1 
C IV (top panel) and O VI (bottom panel) optical depths versus 
H I optical depth for z < 2.4 (filled circles and solid line) and 
z > 2.4 (triangles and dashed line). It can be seen that C IV op- 
tical depth shows mild cosmic evolution when no evolution is 
detected in tovj. The most interesting feature however is that 
the excess of O VI optical depth at log thi ^ —0.5 is more im- 
portant in the high redshift bin. Note that this signal is due to 
the three lines of sight toward quasars with the highest emission 
redshifts. Moreover, if we apply the same procedure as previ- 
ously, we find that, at z < 2.4, the O vi absorption does not 




-1.0 -0.5 0.0 0.5 

lOg(THl) 



Fig. 12. Same as figure 10 for the two redshift bins, z < and > 2.4. 

seem to be related to strong lines when the relation is apparent 
at larger redshift (see Fig 12). As this signal is probably related 
to galactic winds, this may be an indication that galactic winds 
are more prominent at higher redshift. 

5. Conclusion 

We have used a pixel-by-pixel analysis to investigate the pres- 
ence of metals in the inter galactic medium. We have measured 
the median optical depth of C IVA1548 and O VIA1031 absorp- 
tions versus the H I Lyman-a optical depth of the gas at a mean 
redshift z = 2.6 using 19 lines of sight observed at high spectral 
resolution (R = 45000) and high and homogeneous S/N ratio 
(approximately 35 and 70 per pixel over, respectively, the O VI 
and C IV wavelength ranges). Great care has been exercized to 
determine the continuum. 

We have carefully determined the H I optical depth by tak- 
ing into account the information in the Lyman series (usually 
Lyman-/3 and Lyman-7). We have also carefully removed from 
the spectra all wavelength ranges that are polluted by strong 
intervening absorptions from metal line systems, in particular 
Mg II and Fe 11 systems, and from associated systems. 
We find that the gas is enriched in carbon and oxygen for 
Thi > 1- Contrary to previous claims, there is no indi- 
cation that C IV absorption is statistically associated with 
gas of H I Lyman-a optical depth smaller than 1. In ad- 
dition, our observations strongly suggest that the C IV/H I 
ratio decreases with decreasing thj. We observe that for 
Thi ^ 1, log tciv/thi < —3.2 which corresponds to 
log C iv/H l< —3.3 assuming 6hi/&civ = 2 and a Voigt pro- 
file or log C iv/H l< —3.0 assuming 6hi — bcw- This does 
not prevent a small fraction of the low thi gas from being as- 
sociated with strong metal lines (see Fig. 4) indicating that en- 
richment is highly inhomogeneous. 

We detect the presence of O VI for thi > 0.2, consitent 
with a constant ratio log tqvi/thi ^ —2.0 coiTesponding to 
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log O vi/H I ~ —1.75 assuming bm/bovi = 2 and a Voigt 
profile or log O Vi/H i ~ —1.45 assuming 6hi = &ovi • We 
show that for 0.2 < thi < 1, the O VI absorption is associated 
with gas located within ~300 km s~^ from strong H I lines. 
This suggests that the O vi phase is probably part of winds 
flowing away from overdense regions. 

It is therefore not possible to conclude that metals are present 
in the most teneous regions of the IGM (thi < 1 corresponding 
to overdensities of about 1 to 3 at 2; ~3 and 2 respectively) far 
away from overdense regions. 
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